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Abstract 

Meson loop corrections to baryon axial currents are computed in 
the 1/N expansion. It is already known that the one-loop corrections 
are suppressed by a factor 1/N; here it is shown that the two-loop cor- 
rections are suppressed by 1/N 2 . To leading order, these corrections 
are exactly what would be calculated in the constituent quark model. 
Some applications are discussed. 
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1 Introduction 



How are the baryons' properties renormalized by pion loops? This classic 
question gains renewed interest with the advent of each new calculational 
technique. 

Pion loop corrections to baryon properties have been studied using the 
non- linear sigma model with derivative couplings [|IJ. Later, Jenkins and 
Manohar 0, || |J simplified the problem by invoking the heavy baryon ap- 
proximation ||. Using a Lagrangian that included the baryon octet and 
Goldstone boson octet, they found that the one- loop correction to the baryon 
axial current was large — as much as 100% of the tree-level value. However, 
if the baryon decuplet is also included, the total one-loop corrections are 
smaller, on the order of 30% of the tree-level value. That is, the loops involv- 
ing decuplet states tend to cancel the loops involving only octet states. This 
was good news for perturbation theory, but it left unanswered the question, 
"What is the loop expansion parameter?" A seemingly coincidental cancel- 
lation of large corrections did not leave behind any obvious parameter that 
could justify, for example, the belief that the two-loop corrections should be 
any smaller than the one-loop corrections. 

This question can be addressed within the framework of large-N tech- 
niques for baryons || [7|, which have recently been rediscovered and greatly 
expanded || |, [ljj, 0, |]| H (T5J . One of the results is that the baryon 



coupling to the axial current is on the order of the number of colors, 

9A ~ N (1) 



This also means that the baryon-pion coupling is ~ Qa^^I / ~ y N . However, 
the renormalization graph of Fig. 1 (a) gives a contribution of order N 2 
to Qa-i which if taken alone would violate Eq. (|I|) and doom perturbation 
theory (since the one-loop contribution would be much larger than the tree- 
level value). But there is another graph that must not be forgotten, the 
wavefunction renormalization of Fig. 1 (b). When both of these diagrams 
are included, it is found || |ItJ that the leading order behaviour cancels, 



and the total one-loop correction is 0(1), or 1/N times the tree-level value. 
The one- loop corrections are therefore small in the 1/N expansion and chiral 
perturbation theory seems to be valid. 

As encouraging as this result is, it leaves some questions open. If the one- 
loop results are to be fitted to the data and believed, one should show that 
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the two-loop contribution is small compared to the one-loop result. This is 
not obvious since, for example, the diagram of Fig. 2(a) is of order N 3 times 
the one-loop correction. However, once again there are several diagrams to 
be added together. This paper shows that when all of the two-loop diagrams 
are taken into account, the largest terms cancel, and the result is of order 
1/N times the one-loop contribution. Evidently, the pion loop expansion 
parameter turns out to be 1/N. Although this has been suspected before 
[§, it has not been previously demonstrated to two loops. 



One result of this analysis is a demonstration that when the pion-baryon 
vertex is taken to leading order in 1/N, the chiral loop corrections are the 
same as in the chiral quark model. This is surprising because in that model 
[ PH , |3|] the quarks are considered to be weakly bound, so that the pions 
interact with only one quark at a time, and the other quarks are spectators. 
In the present 1/N analysis, however, the quarks are so tightly bound that 
the pion interacts with the baryon as a whole, and there are no spectators. 
Nevertheless, when all of the diagrams are taken into account, many cross- 
terms cancel, and the 1/N result is just what would have been calculated 
in the chiral quark model. It had already been noted that the constituent 
quark model fits the data as well as the usual baryon-pion theory ||; the 
1/N expansion sheds some light on why this is so. 

The next section presents the two-loop calculation; this is followed by a 
discussion of possible applications of this formalism. 



2 Two-Loop Corrections 

What are the meson loop corrections to the baryon axial current? The spatial 
components of the axial current are written M, [10], |TTJ 

<&\q<y i <y B T a q\B> lTee =X% B 

where T a is a generator of the flavor group. (The magnitude of X ta is what 
was loosely called gA in the introduction.) The axial current can be expressed 
in the 1/N expansion using perturbative baryon states \B > [0 (to be 
contrasted with physical states \B >): 

\B >= B riai --- TNaN V o t 10 > 
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Here a\ a creates a quark with spin a and isospin r. The quarks are totally 
antisymmetric with respect to color; their color indices are suppressed, and 
the operators a and a) are treated as bosonic, rather than fermionic [|T^]. 
Then the axial current can be written in a 1/N expansion 

X% B = g < B'\G ia \B > +j^< B'\H ia \B > + ... (2) 
where g and h are constants of order 1 and G ta and H ta are the operators 



12, 13, 15 



Gia f rpa i 

— a ra ± rs a aj3 a s /3 

H ia = (aJ Q ,T r a s a sa )(aJ /3 cr^ 7 a 47 ) 

The mesons are coupled derivatively to the baryon axial current. Some 
of the Feynman rules for the pion-baryon interactions are given in Fig. 3. 
The baryons are treated within the heavy fermion approximation |5], ^|, [J, 
and the calculations are performed in the baryon's rest frame. The meson 
propagator uses the mass matrix m 2 b , a diagonal matrix that gives the masses 
of the pions, kaons, and eta under flavor symmetry breaking. For the N 
power counting, it is important to keep in mind that the pion decay constant 
focVN. 

Now look at the vertex renormalization. The momentum integral for Fig. 
1(a) is 

d 4 p 1 |p 2 



Z oh = 'l I 777Z\4 ^2 



(2vr) 4 pi p 2 - m\ h 
For Fig. 2(a) and Fig. 2(c), the integral is 



J? 



d 4 p d*q 1 f 1 1 \ |p 2 |q 2 

(2tt) 4 (2tt) 4 pi \p + q Q q J p 2 - m 2 aa , q 2 - m 2 w 



The inner loop includes a counter-term. If this term (the 1/go appearing 
above) is not included, the internal baryon acquires an additional mass, which 
must then be transformed away by the heavy baryon transformation. It is 
easier to simply include the counter-term explicitly. The mass differences 
between the various baryons are proportional to 1/N and/or to the flavor 
symmetry breaking, and will be ignored. 
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The integrals for Fig. 2(b),(d),(e), and (f) are (respectively) 

rraa'bb' / 

J 2 J 



d 4 p d A q 1 1 |p 2 |q 2 

WY A WT 4 Pi (po + Qo) 2 P 2 - < a , q 2 ~ m 2 w 



w , i_ i i |p 2 |q 2 



(2tt) 4 (2vr) 4 p 2 , po + go <?o P 2 - 9 2 - ™ 2 fe ' 
_ / #P dS 1 1 1 |p 2 |g 2 



(2tt) 4 (2tt) 4 po (po + <?o) 2 ?oP 2 - ^L' <? 2 ~ m w>' 

/ _d^_d*q_ J_ 1 1 |P 2 _ £ 66 'a,/ 

* (2tt) 4 (2vr) 4 po po + q q 2 p 2 - m\ a , g 2 - m 2 bb , ! 
The vertex renormalization to two loops can then be written: 

V£ B = \X ia + Lxbb' x jb' x ia x jb + J_jbb>cc> X jb x ia x kc x kc> X jb' 

V / / 

j^^j-bb'cc' X jb x kc x ia x kc' x jb' _|_ sj-bb' cc' x jb x kc x kc' x ia x jb' 

I - cc' x jb x ia -\rkc x jb' -\rkc' _|_ ^ j^bb'cc' -\rjb -\rkc x ia x jb' X kc' 

/4 1 f 2 

+ 1 K <x>bb> x jb x kc X jb> X ia x kc> j ( 3 ) 



/ 4 



B'B 



The operators X M are treated as matrices with baryon indices; intermediate 
baryon states are summed over. 

The baryon wavefunction renormalization constant can be computed from 
the diagrams of Fig. 1(b) and Fig. 4: 

(Z- 1 ) bib = (l + ^f bb 'X jb X jb ' + j Z (2J* bfccf + jWcc^ x jb x kc X kc' xj b' 
+71 W cc ' + ^2^' + bb ')X jb X kc X jb 'X kc ' j (4) 
Finally, the renormalized axial current is 



i i 



< B'\qf l5 T a q\B >= Zj V ia Z% (5) 

B'B 
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When Eqs. (|3|) and ([|) are substituted into Eq. (|) and the result is 
multiplied out to order l// 4 (i.e. to two loops), the terms do not simplify in 
any obvious way. Some identities among the integrals must be used, 



aa'bb' 



Y*aa'bb' 

/c 2 



aa'bb' 



jQaa'bb 1 _|_ sji 



aa'bb' 



j^aa'bb' _|_ j^bb'aa' -j-aa' 



bb'aa' 



raa' -rbb' 



as well as the identity 



j^aa'bb' ^y^ia j^ia' X^X^ 



Then after a few pages of algebra, Eq. (|5p becomes 

< B'\qYl 5 T a q\B >= X ia + -Ll bb '[X jb , [X ia ,X jb '}} 

2f 2 



I 1 trbb'cc' 



1 



X Jb , [[X kc , [X ia ,X kc '}],X J 



ia vkc'l 



+ 2 



[X jb ,X% [X kc , [X jb ',X kc '}} 



+ 



X 1 



4/' 



■K, 



bb'cc' 



X j \ [X kc , \x 3h \x kc ' 
[X jb ,X kc ],[X ia , [x jb ',x kc> 



(6) 



This is the main result of the paper. From here it is possible to show that 
the one-loop corrections to the axial current are suppressed by 0(1/ N) times 
the tree-level, and the two-loop contribution is suppressed by 0(1/N 2 ). 
For example, suppose we take X ia to leading order in 1/N, 



X ia = g< B'\G ia \B > 



(7) 



where G m = a\ a T^ s a l a pa s ^. Since the operator G m has one a and one a\ it 
can count the number of quarks in the baryon once, and can therefore be of 
order N . No accidental cancellations occur, so the axial current is O(N) to 
tree level. Now using Eq. ([?[), the one- loop correction can be read off from 
Eq. ©: 



jjf bb ' < B'\\G 3b , [G ia , G jb ']]\B > 



(8) 
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Each G ia has one a and one a', but each commutator eliminates an a-a) pair 
due to the identity [al a V^a b p, al^W^a d s} = aLjV, W] b a ^abp- The resultant 
operator in Eq. (§) has only one a and one a', so the matrix element is at 
most of order N. Since 1/f 2 ~ 1/iV, the total one-loop correction is 0(1), 
or 1/JV times the tree-level value. 

The quadruple commutators, which give the two-loop corrections, also 
eliminate all but one a-a^ pair. Therefore these commutators are also of 
O(N), and when they are multiplied by a coefficient of l// 4 , the result is 
0(1/N). That is, the two-loop correction is 1/N 2 times the tree-level value. 

This result has the following interpretation: since the loop corrections 
contain only one a and one a', the pion vertices and the current operator 
all act on the same quark; the vertex and wavefunction renormalization are 
carried out on each quark individually. 

The reason that this happens in the one-loop case can easily be demon- 
strated graphically. In Fig. 5(a), the ends of the pion line are attached to 
different quark lines than the one on which the current acts. In such a case, 
one of the pion vertices can be commuted past the current operator, so that 
the pion is emitted and absorbed before the current operator acts (see Fig. 
5(b)). However, this diagram has already been taken into account by the 
wavefunction renormalization, and so does not contribute. Similarly, in the 
graph of Fig. 5(c), the pion line begins on a different quark line from the one 
on which the current acts, but ends on the same line. In this case, the first 
vertex can be commuted past the current operator, resulting in Fig. 5(d). 
Again, this diagram is cancelled by the wavefunction renormalization. So 
we reach the conclusion stated above: the only type of graph that needs to 
be considered has both pion vertices acting on the same quark line as the 
current operator (Fig. 5(e)). 



This is exactly what is assumed to be true in the chiral quark model |L6 
as developed in Ref. |J. In that model the leading vertex and wavefunc- 
tion renormalizations act on one quark at a time, and the quark propagator 
is that of a (fairly heavy) constituent quark. One problem with the chiral 
quark model, however, is that it is difficult to understand why a free quark 
propagator can be used inside the proton. Why should not the bound state 
wavefunction be taken into account? The large-N approach offers a solu- 
tion: the quark inherits the heavy fermion propagator from the baryon as a 
whole but manages to decouple from the other quarks thanks to the cancel- 
lations that occur in the commutator structure of the loop corrections. The 
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constituent quark emerges from a tangle of meson loops. 

There are two points that complicate the above picture somewhat. The 
first is that we have left out some diagrams. Figs. 6 (a) - (e) also contribute 



they are most easily calculated using an effective Lagrangian [[L2], [13 



[14]]. It turns out that these diagrams follow the same pattern as above: l- 
loop diagrams are suppressed by factors of (1/N) 1 , and the result is just what 
would have been expected from the chiral quark model. One difference of 
the these diagrams, however, is that they are not necessarily suppressed by 
powers of the coupling constant g of Eq.(|7|). For example, Fig. 6(a) and Fig. 
6(b) are both proportional to g (rather than g 3 or g 5 ). This point does not 
affect the present discussion, but is important in the next section. 

The second complication is that Eq. (|7|) is only an approximation to 
the axial vertex. When the vertex is expanded to the next order in 1/N 
(as suggested by Eq. (0) ), a new operator H ia is introduced. H m acts 
on two quarks at a time, so the simplest constituent quark picture receives 
corrections^. The axial current is of O(N), and the leading corrections are 
of 0(1). To evaluate all the corrections to 0(1), the operator H ia must 
be included at tree level, as well as in the double commutator of Eq. @. 
(One might have thought that the loop corrections involving the operator 
H m would be smaller in 1/N than H m itself, but this turns out not to be 
the case.) Similarly, in order to calculate loop corrections to 0(1/N), X ta 
must be expanded to 0{1/N) and then inserted into Eq. (H). Nevertheless, 
the identities of Ref. |y| can be used to show that the one- and two-loop 
corrections are still of (9(1) and 0(1/N) respectively. In actual practice, Ref. 
[ |l3|j has estimated h ph —0.1, and so the operator H ia will be ignored below. 



3 Discussion 

How does all this formalism apply to the real world? The double commu- 
tator in Eq. (P), which gives the one-loop correction, is of the order of the 
number of flavors Np, and the quadruple commutators giving the the two- 
loop corrections are of order Np. The momentum integrals should be cut off 
at the chiral symmetry-breaking scale A. Let the pion decay constant / be 

3 Actually such operators appear in the chiral quark model also; rather than being 
suppressed by 1/N, though, they are suppressed |l6| by a power of the wavefunction at 
the origin divided by the constituent quark mass, |V>(0)| 2 / 3 /m c . 
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factorized to show clearly its N-dependence: 

f = VNf 

where / is 0(1). As mentioned previously, some diagrams of Fig. 6 are not 
suppressed by powers of the axial coupling constants (g and h of Eq. (0)). 
Therefore the chiral loop expansion parameter is 

N F A 2 

167T 2 / 2 

This parameter is not small in any estimation |L7|. However, one can adopt 
the following approach: start with the bare coupling g (or for example h), 
assume that it can be renormalized to all orders in the flavor symmetric limit 
(m v = itlk = vn^ « 0), resulting in the renormalized constant g^. This new 
constant gn is to be used in computations, and the effects of virtual pions 
can be computed loop-by-loop, keeping only those terms that violate SU (N F ) 
symmetry. In this case all terms involving A 2 /16ir 2 f 2 are to be thrown away, 
since their effects have already been included in g^. The new loop expansion 
parameter then becomes 

Nf m\ A 2 
~N~ 16tt 2 / 2 ° g rr& 

This procedure is equivalent to using dimensional regularization for all the 
integrals of Eq. (|5]) . 

Such a program has already been carried out by Ref. ||. The one- 
loop corrections to the baryon axial current were computed in the chiral 
quark model using dimensional regularization. This is exactly equivalent to 
a leading order 1/N calculation. The model fits the data well; a best fit is 
obtained for gn = 0.56. 

So far we have examined the corrections to the octet axial currents. 
What about the singlet current, the "spin content" of the baryon? When 
computing the renormalization of the singlet current, fewer diagrams ap- 
pear than for the octet current: Figs. 6 (a)-(c) do not exist for the sin- 
glet current. Therefore the one- loop contribution is suppressed by a factor 
9r(Nf/N)(A 2 /167T 2 / 2 ) compared to the tree value. The two-loop diagrams 
are suppressed by a factor of g 2 R (N F /N)(k 2 /l§ii 2 f 2 ) (Figs. 2 and 4) or 
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(N F /N)(m 2 K /167i 2 f 2 )(\ogA 2 /m 2 K ) (Figs. 6(d) and (e)) compared to the one- 
loop diagrams. Therefore, the loop expansion parameter e for the singlet 
current is 

( 2 N F A 2 N F m\ , A 2 \ 
e = max g R — — -g- , — -——5- log ■ 



iV 167T 2 / 2 ' AT 16tt 2 / 2 
If e is small enough, we do not have to go through the extra step of first do- 
ing the flavor-symmetric renormalization and then returning to the integrals 
using dimensional regularization. Cutoff regularization can be used from the 
outset. 

Using this approach, and the leading 1/N baryon-pion vertex, the spin 
content of the proton turns out to be 




32 A 2 22 m 2 K , (2A) 2 

v log 



9 16vr 2 f 2 9 167r 2 f 2 m 



K 



(9) 



(Here I used m 2 = and m 2 = jW^.) This calculation assumes that A 
is a physically meaningful cutoff (it is the point where chiral symmetry is 
restored). This assumption differs from other chiral perturbation calculations 
0, 0, £|, where the quadratically divergent pieces are discarded and 
only the logarithmically divergent pieces are kept. However, retention of 
the quadratically divergent pieces is well within the spirit of the effective 
Lagrangian Jl6| . Unfortunately, since neither g nor A is known, this equation 



has no predictive power. It is comforting, however, that a reasonable choice of 
the parameters gives a reasonable result. For example, for g = 1 and A = 1 
GeV, Eq. (^j yields a spin content of 0.57, which is within 0(e 2 ) of the 
experimental value |18| of 0.27 ±0.11 . In this case the expansion parameter 



e is rather large, e ~ 0.75, so the effects of chiral loops are estimated to be 
very important. 
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Figure Captions 



Fig. 1 The one loop vertex renormalization (a) and wavefunction renor- 
malization (b). The "x" represents the axial current operator. 

Fig. 2 Two loop vertex renormalization diagrams. 

Fig. 3 Some of the Feynman rules for baryon-meson interactions. 

Fig. 4 Two loop wavefunction renormalization diagrams. 

Fig. 5 The renormalization of the baryon axial current, viewed at the 
quark level. Diagram (a) = diagram (b), and (c) = (d). Only (e) ends up 
contributing to the overall renormalization. 

Fig. 6 Additional one and two loop diagrams that contribute to the 
renormalization of the axial current. 
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